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Abstract 

This  paper  introduces  a  general  and  powerful  framework  for  the  analysis  of  uncertain 
systems,  encompassing  linear  fractional  transformations,  the  behavioral  approach  for  sys¬ 
tem  theory  and  the  integral  quadratic  constraint  formulation.  In  this  approach,  a  system 
is  defined  by  implicit  equations,  and  the  central  analysis  question  is  to  test  for  solutions 
of  these  equations.  In  Part  I,  the  general  properties  of  this  formulation  are  developed, 
and  computable  necessary  and  sufficient  conditions  are  derived  for  a  robust  performance 
problem  posed  in  this  framework. 


1  Introduction 

In  the  predominant  viewpoint  in  systems  and  control  theory,  a  system  is  an  input-output  (I/O) 
entity,  where  the  variables  are  clearly  separated  in  two  groups,  and  a  cause-effect  relationship 
is  established  between  them.  This  approach  entails  a  “signal  flow”  conception,  adequate  for 
systems  which  are  deliberately  built  to  match  the  I/O  philosophy,  such  as  computers  and 
amplifiers. 

For  many  other  physical  systems  this  point  of  view  often  appears  artificial;  as  an  example, 

a  mass  or  energy  balance  equation  in  a  chemical  process  is  more  naturally  thought  of  as  an 

equation  or  constraint  between  variables  than  as  a  cause-effect  law.  While  this  observation 
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will  appear  as  no  surprise  to  an  engineer  performing  modeling  of  such  a  system,  it  is  only 
recently  that  its  theoretical  implications  have  been  extensively  considered. 

In  a  research  program  best  summarized  by  the  survey  paper  [19],  Willems  has  advocated 
an  approach  to  system  theory  where  the  central  concept  is  the  behavior,  a  set  of  allowed  signal 
trajectories,  and  no  input-output  partition  is  a  priori  established  between  the  variables.  The 
corresponding  theory  of  finite  dimensional  linear  systems  has  been  extensively  developed  [19] . 

This  paper  adopts  the  same  philosophy  and  contends  that  this  point  of  view  is  even  more 
natural  for  systems  involving  uncertainty.  If  the  relationship  between  the  variables  is  not 
precisely  known,  the  cause-effect  point  of  view  is  itself  suspect:  it  is  more  natural  to  think  of 
an  uncertain  implicit  relationship  between  variables. 

It  is  noteworthy  that  partial  versions  of  this  viewpoint  have  been  present  in  early  work 
leading  to  current  robust  control  theory.  In  the  foundational  paper  [22],  Zames  states  the 
basic  stability  theorems  describing  systems  as  relations,  motivated  by  some  nonlinear  systems 
which  do  not  fit  the  I/O  concept;  the  same  ideas  are  present  in  Safonov  [16]  in  the  early  years 
of  robust  control.  A  further  example  which  has  led  to  powerful  analysis  techniques  is  the 
Integral  Quadratic  Constraint  (IQC)  formulation  of  Yakubovich  [21]  and  Megretski  [9,  10], 
where  a  component  is  described  by  constraints  between  the  signals  involved. 

It  might  be  argued  that  as  long  as  the  I/O  framework  involves  no  mathematical  loss  of 
generality,  an  argument  of  symmetry  or  esthetics  is  not  sufficient  to  abandon  an  approach 
which  is  widespread.  This  paper  shows  evidence,  however,  that  with  respect  to  robustness 
analysis,  an  implicit  approach  strictly  enhances  the  range  of  applications  of  existing  theory 
for  the  I/O  setting.  The  main  extension  which  is  provided  is  the  ability  to  analyze  over¬ 
constrained  systems  in  a  unified  framework;  these  arise  when  superimposing  an  nncertain 
model  and  a  number  of  constraints  relevant  to  the  analysis  problem  under  consideration. 

In  Part  I  of  this  paper  we  propose  a  theoretical  framework  which  encompasses  the  behav¬ 
ioral  approach  for  system  theory,  the  Linear  Fractional  Transformation  (LFT)  paradigm  for 
uncertainty  descriptions,  and  the  IQC  formulation.  Section  3  introduces  the  framework  and 
shows  how  it  allows  for  the  formulation  of  a  general  robust  performance  analysis  problem,  in 
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terms  of  finding  solutions  to  uncertain  equations. 

This  problem  leads  to  a  notion  of  robust  stability  for  implicit  systems,  presented  in  Section 
4,  which  naturally  extends  the  existing  I/O  theory.  The  general  analysis  problem  is  reduced 
to  a  canonical  case  to  be  considered  in  the  rest  of  the  paper. 

Section  5  contains  an  analysis  test  for  robust  stability  of  an  implicit  system  when  the 
nominal  equations  are  linear,  time  invariant,  and  the  uncertainty  is  allowed  to  vary  in  the 
class  of  arbitrary  norm  bounded  operators.  This  condition  is  a  convex  feasibility  test  on 
a  constant  scaling,  which  extends  the  scaled  small-gain  conditions  for  robust  stability,  and 
recent  results  [17,  10]  on  the  necessity  of  these  conditions  for  the  standard  I/O  setting.  The 
extension  also  includes  61  operator  blocks  in  the  uncertainty  description. 

Preliminary  versions  of  these  results  were  presented  in  the  conference  papers  [14,  15].  All 
proofs  are  collected  in  the  Appendix. 


2  Notation 


In  this  paper  we  will  consider  vector  spaces  of  signals  V  C  (F")^,  where  F  is  the  field  R  or 
C,  and  the  time  index  T  can  denote  continuous  time  R  or  R"*",  or  discrete  time  Z  or  Z+;  for 
concreteness,  the  results  will  be  presented  in  discrete  time,  but  most  of  the  theory  extends  to 
continuous  time  in  a  straightforward  way.  For  T  £  T,  define  the  truncation  operator  as: 


P’’  :  v{t)-^ 


v(t)  if  t  <  T 
0  if  f>T 


(1) 


Signal  norms  can  be  introduced,  choosing  V  as  a  Banach  space  of  signals.  Following 
[22,  5],  we  introduce  extended  signal  spaces  to  include  signals  which  “blow  up”  at  infinite 
time,  by  defining  (for  positive  time  axis  T  =  Z+)  —  {v  E  (F”)^  :  P^{v)  G  V  VT  G  T}. 

In  this  paper  we  will  consider  the  2-norm  for  signals;  we  will  indicate  later  which  parts  can 
be  extended  to  other  signal  norms,  will  denote  the  Hilbert  space  of  square  summable, 
F"-valued  sequences  over  Z  or  Z+;  ||u||^  =  XlteT  ^^.e  corresponding  extended 

space  over  Z"*".  The  spatial  dimension  n  will  be  dropped  when  clear  from  context. 

The  class  of  linear,  bounded  operators  G  :  i®  denoted  or  simply  C{l2).  We 
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are  also  interested  in  operators  which  can  be  extended  to  l2e',  in  this  respect  it  is  convenient  to 
consider  the  class  of  causal m&ps,  which  verify  P^GP^  —  P'^G  for  all  T.  Equivalently,  a  causal 
linear  map  can  be  characterized  by  having  a  lower  triangular  infinite  matrix  representation. 

A  causal  linear  map  over  /2e  has  finite  gain  if  there  exists  7  <  oo  such  that  ||P^Gu||  < 
7  ||P^n||  for  all  V  G  ke,  T  eJ.  Equivalently,  the  restriction  of  G  to  I2  is  a  bounded  operator, 
of  norm  ||G||  =  supyg2+(^  (^^))  where  G’^  is  the  matrix  representation  of  P^GP^ , 


G'^ 


Goo  0  0 

G'lo  Gii 


(2) 


Gj’o  Gxi  ■  ■  ■  G'j't 

and  a  denotes  maximum  singular  value.  We  will  also  use  the  minimum  singular  value  of  a 
matrix  ct(G')  =  min{||G^||  :  ||^||  =  1}. 

The  set  of  all  causal,  finite  gain  operators  G  :  denoted  simply 

Cc{he)-  The  unit  delay  operator  is  denoted  by  A.  A  map  G  is  time  invariant  if  G\  =  AG. 


3  The  Implicit  Framework  for  Analysis 


3.1  Implicitly  Defined  Systems 


This  paper  deals  with  implicit  characterizations  of  systems.  Loosely  speaking,  this  means  that 
the  laws  governing  a  system  and  the  constraints  imposed  on  a  problem  under  consideration 
are  all  expressed  as  equations  on  a  specified  set  of  variables.  A  formal  definition  follows. 

Definition  1  An  implicit  system  (W,E,G)  is  defined  by  two  vector  spaces,  the  variable 
space  W  and  the  equation  space  E,  and  an  equation  operator  G  :  W^E.  The  behavior 
of  the  implicit  system  is  the  set  B  =  Ker{G)  =  {w  £  W  :  Gw  =  0}.  The  system  is  called 
linear  if  G  is  a  linear  map. 

The  definition  above  is  closely  related  to  the  behavioral  approach  to  system  theory,  introduced 
by  Willems  [19].  In  this  type  of  formulation,  all  variables  in  a  system  are  a  priori  on  an  equal 
footing,  without  a  distinction  between  inputs  and  outputs.  The  system  laws  are  constraints 
in  the  possible  values  of  these  variables,  which  define  a  set:  the  behavior. 
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In  our  case  the  implicit  equations  defined  by  G  play  a  central  role,  not  captured  entirely 
by  the  behavior,  since  as  we  shall  see  equation  error  will  be  added  for  the  analysis.  Intercon¬ 
nections  of  subsystems  is  reduced  to  superimposing  equations. 

These  descriptions  arise  naturally  when  modeling  physical  systems  from  first  principles, 
where  physical  laws  are  more  naturally  thought  of  as  equations  between  variables  than  as 
“signal-flow”  maps.  For  instance,  a  resistor  in  a  circuit  is  naturally  modeled  by  the  equation 
v  —  Ri  =  0,  and  there  is  no  need  to  specify  a  cause-effect  relationship  between  the  two  variables 
V  and  i;  this  distinction  is  artificial  and  not  available  a  priori.  For  further  discussion  of  the 
features  of  this  modeling  paradigm,  see  [19,  2]. 

An  important  special  case  of  Definition  1  is  the  class  of  dynamical impYiciX  systems,  where 
the  sets  W  and  E  are  vector-valued  signal  spaces.  As  an  example,  if  R{^)  is  a  polynomial 
matrix,  the  differential  equations  R{j^)w  =  0  define  an  implicit  system,  where  G  is  the  differ¬ 
ential  operator  R{-^),  and  W,  E  can  be  chosen  as  spaces  of  smooth  functions  (or  alternatively 
distribution  spaces). 

The  choice  of  the  defining  elements  W,  E,  G  is  essentially  determined  by  the  type  of 
analysis  to  be  performed.  In  this  paper  we  will  consider  two  possible  settings  which  are 
standard  in  robust  control: 

1.  To  formulate  quantitative  performance  specifications,  a  signal  normed  space  is  required; 
we  use  the  I2  space  W  =  ll,  E  =  I2,  with  G  a  bounded  mapping  between  them.  This 
paper  deals  with  linear  equations  G  G  E{h)- 

2.  In  some  cases  in  stability  theory  it  cannot  be  assumed  a  priori  that  signals  have  finite 
norm,  and  the  extended  spaces  W  —  ll^,  E  =  /fe  ^re  used,  together  with  G  G  Cc{he)- 

These  two  cases  will  be  followed  in  parallel  throughout  this  paper. 

3.2  Uncertainty  and  LFTs 

We  now  incorporate  into  the  implicit  paradigm  deterministic  descriptions  of  uncertainty  in  the 
style  of  robust  control.  The  equation  map  G  is  replaced  by  a  parameterized  map  G(A),  where 


5 


A  is  an  uncertainty  operator.  A  general  class  of  uncertainty  descriptions  can  be  parameterized 
by  the  class  of  “spatially”  structured  perturbations,  of  the  form 

A  =  diag  Silr^,  A^^i, . . A^+j’]  (3) 

The  diag  notation  implies  that  Az  =  [SiZi,...,6lZl,...Alj^pZl+f\',  where  the  vector  z  is 
broken  into  spatial  components  Zi  of  appropriate  dimensions.  The  blocks  Ap+j  can  be 
used  to  describe  real  parameters,  or  dynamic  (linear  time  invariant  (LTI),  linear  time  vary¬ 
ing  (LTV)  or  nonlinear)  perturbations;  the  square  blocks  correspond  to  an  uncertainty 
perturbation  which  is  repeated  in  the  spatial  dimension;  the  possibly  non-square  full  blocks 
Ai+j  are  unrestricted  maps.  In  each  case,  there  is  a  restricted  class  A  of  allowed  pertur¬ 
bations.  In  this  paper  we  restrict  ourselves  to  linear  uncertainty  blocks,  i.e  A  C  ^(/j)  or 
A  C  £c(^2e)-  The  uncertainty  will  be  normalized  to  a  unit  ball  Ba  in  the  operator  norm, 
Ba  -  {A  G  A  :  ||A||  <  1}. 

For  the  parameterization  <j(A)  of  implicit  uncertain  systems  we  will  adopt  the  linear 
fractional  transformation  (LFT)  paradigm  (see  [11])  depicted  in  Figure  1,  which  provides  rich 
uncertainty  descriptions.  Implicit  LFT  representations  were  first  considered  in  [3]. 


_ m _ 

-  A  B  *- 

—  2: 

0^ 

^  C  D  ^ 

—  w 

Figure  1:  Implicit  LFT  system 


In  Figure  1,  M  =  ^  or 

o  u 

portant  special  case  is  M  time-invariant.  The  uncertainty  A  has  the  structure  (3). 

A  remark  regarding  Figure  1  is  that  it  contains  remnants  of  the  “signal-flow”  approach, 
since  the  parameter  A  is  depicted  as  an  input-output  map.  This  is  done  to  highlight  the 
connection  with  the  standard  LFT  paradigm,  but  from  an  implicit  point  of  view  the  system 
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in  Figure  1  is  simply  characterized  by  the  equations 

<p{A,M)  ^  =0,  where  (p(A,M)=  ^  ^  (4) 

This  description  is  “internal”  since  the  signals  z  produced  by  the  uncertainty  operators  are 
included  in  the  variable  space  W.  The  parameterization  G[A)  is  therefore  affine  in  the 
parameter  A.  This  simple  form  allows  the  representation  of  a  rich  variety  of  uncertain  systems; 
in  fact,  all  the  complexity  is  captured  by  the  structure  of  A.  A  standard  input-output  LFT 
uncertain  system  can  be  easily  converted  to  this  implicit  form  (see  section  3.4  below). 

3.3  Integral  Quadratic  Constraints 

We  will  now  present  a  feature  of  the  implicit  analysis  framework  which  does  not  have  a 
counterpart  in  the  I/O  setting:  it  allows  for  the  representation  of  additional  constraints  in 
the  signals  of  a  robustness  analysis  problem. 

In  particular,  the  implicit  formulation  over  I2  permits  the  representation  of  Integral 
Quadratic  Constraints  (IQCs),  which  have  been  proposed  by  Yakubovich  [21]  and  Megret- 
ski  [9]  as  the  basis  of  an  alternative  paradigm  for  robust  control.  IQCs  are  time-invariant 
quadratic  forms  in  signal  space,  which  in  discrete  time  have  the  form 

{Jiz,z)=  r  z{e?'^yn{e^'^)z{e^'^)dLo  <0  (5) 

J  —  TT 

where  z(e^'^)  is  the  Fourier  transform  of  an  I2  signal,  II  =  11*  is  an  LTI  operator,  assumed 
bounded  on  I2  (i.e  n(e-^“')  is  in  Loo)- 

IQCs  can  be  used  to  provide  deterministic  descriptions  of  an  uncertain  component,  by 
defining  a  set  of  signals  which  captures  the  known  information  about  the  component  (as  in 
[9]).  As  will  be  explored  in  Part  II  of  this  paper,  it  may  be  of  interest  to  describe  properties  of 
a  disturbance,  by  constraining  it  in  terms  of  IQCs;  of  siich  nature  are  “whiteness”  constraints 
as  in  [12],  which  describe  spectral  properties  of  a  disturbance. 

It  is  now  shown  that  these  general  IQCs  can  be  captured  by  an  uncertain  implicit  system. 
First  choose  P  LTI  (e.g.  P  =  kl)  such  that  P* P  —  11  >  0.  A  spectral  factorization  gives 
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p(giw)*p(gjw)_II(gjw)  _  Q(^e^‘^)*Q(e^'^),  where  Q  can  be  chosen  in  Hoo-  So  11  =  P*P  —  Q*Q, 
which  reduces  (5)  to  II-P2II2  <  We  now  introduce  the  following  lemmas: 

Lemma  1  Let  z  G  /”,?;  €  The  following  are  equivalent: 

(i)  11^112  -  Ikll2  >  0 

(n)  3A  G  JC{12,1?),  ||A||  <  1  :  Av  =  z 

Lemma  2  Let  z^v  E  l^-  The  following  are  equivalent: 


(i)  /  u(e'")u(e-'")*  —  z(e^‘^)z(e^‘^)*  du  >  0 

J  —Tt 

(8) 

(m)  V77  G  CW\7]*v\\^  >  \\rfz\\^ 

(9) 

(in)  36  G  £(/2),  ll^ll  <  1  :  hL^v  =  z 

(10) 

(6) 

(7) 


From  the  above  discussion  and  Lemma  1,  for  each  2;  satisfying  (5)  there  exists  an  operator 
Ac,  ||Ac||  <  1  such  that  Pz  =  AcQz.  So  the  set  of  2:  G  ^2  satisfying  (5)  can  be  described  as 
the  union  over  Ac,  ||Ac||  <  1  of  the  behavior  of  the  uncertain  implicit  system 

{P-AcQ)w  =  Q  (11) 


Remarks: 

•  Although  for  each  Ac  the  set  Ker{P  —  Ac<5)  is  a  linear  subspace,  the  union  of  the 
parameterized  behaviors  describes  a  more  complicated  set  given  by  (5). 

•  The  set  ||Ac||  <  1  includes  arbitrary  time- varying  non-causa/ operators.  In  this  respect, 
implicit  systems  obtained  from  this  procedure  are  a  priori  considered  in  the  I2  setting. 

•  A  finite  number  of  IQCs  can  be  given  a  representation  (11),  where  Ac  is  now  a  struc¬ 
tured  uncertainty  operator. 
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•  Analogously,  the  set  U||Ac||<i  Ker{P  —  del  Q)  for  scalar  6c  can  be  shown  by  means  of 
Lemma  2  to  correspond  to  the  matrix-valued  constraint 


r  P(e-'“)2(e-'“)2(e-'“)*P(e^“)*  -  du  <  0  (12) 

J  —  TV 


Matrix  valued  constraints  appear  naturally  in  the  context  of  characterizing  multivariable 
white  noise  disturbances,  as  shown  in  Part  II. 


3.4  Robust  performance  analysis  in  the  implicit  setting 

To  illustrate  how  implicit  descriptions  might  be  used  for  robust  performance  analysis,  consider 
the  uncertain  I/O  system  of  Figure  2,  where  it  is  known  that  the  (disturbance)  input  d  satisfies 
a  finite  number  of  restrictions  in  terms  of  IQCs  as  in  (5).  We  want  to  determine  whether  there 
exist  signals  d  in  the  allowed  class,  and  perturbations  such  that  the  system  I2  gain  is  1  or 
larger.  This  last  requirement  is  captured  by  the  extra  “performance  IQC”  ||d||^  -  \\y\f  <  0. 


1 

Pii  Hu 

y  - — - 

H21  H22 

Figure  2:  Input/Output  LFT  system 


The  implicit  equations  for  the  system,  the  IQCs  on  d  and  the  performance  IQC  are 
captured  respectively  by  (13), (14)  and  (15),  where  Ac  and  Ap  are  norm  bounded  operators 
(Ac  is  in  general  structured). 


■/-A„Pn 

0 

-A„Px2  ■ 

P21 

-I 

H22 

y 

d 


=  0 


(P- Ac(3)d  =  0 


(13) 
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(14) 

(15) 


The  superposition  of  (13),  (14)  and  (15)  gives  an  implicit  description  for  the  robust  per¬ 
formance  analysis  problem,  which  essentially  reduces  to  the  question: 

Q:  “Does  there  exist  a  perturbation  (Au,Ac,Aj>)  such  that  (13-15)  have  non¬ 
trivial  I2  solutions?”. 

In  Q  the  input/output  partition  has  been  eliminated  from  the  problem,  and  the  analysis 
is  posed  in  terms  of  equations  and  solutions,  rather  than  maps  and  gains;  this  allows  for  a 
natural  incorporation  of  the  constraints  (14).  Questions  of  this  type  are  analyzed  in  the  rest 
of  this  paper. 


4  Stability  in  Implicit  Systems 

We  begin  by  reviewing  the  concept  of  stability  in  standard  system  theory.  Referring  to  the 
M-N  feedback  interconnection  of  Figure  3  (a),  stability  can  be  given  two  interpretations. 


Figure  3:  Stability  in  a  standard  feedback  interconnection 


In  the  first  place,  from  the  point  of  view  of  dynamical  systems,  stability  ensures  that 
solutions  do  not  exist  to  the  loop  equations  where  the  signals  are  unbounded  (e.g.  in  heXh)- 
This  notion  is  usually  considered  for  causal  systems  M,  N. 

Secondly,  from  an  operator  theoretic  point  of  view,  stability  is  sometimes  interpreted  as 
bounded  sensitivity  to  small  errors  or  disturbances  (ei,  62  as  in  Figure  3  (b))  injected  at  the 

r  / 

interconnection.  In  other  words,  the  map  ^  ^  between  Ci,  and  signals  w,z  is  a 

bounded  operator  in  some  normed  space.  This  notion  can  be  stated  for  non-causal  systems 
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M  and  N,  provided  all  signals  are  a  priori  constrained  to  the  normed  space  (e.g.  I2)  where 
the  operators  M  and  N  are  well  defined. 

Although  these  two  versions  are  equivalent  in  many  special  cases  of  causal  linear  systems, 
we  will  find  it  useful  to  distinguish  the  two  for  the  extension  to  the  implicit  framework,  since 
we  are  led  naturally  to  include  non-causal  perturbations  as  explained  in  Section  3.3.  We  will 
term  the  first  notion  “stability”  and  the  second  “/2-stability”.  All  the  material  in  this  section 
can  be  extended  to  other  signal  norms  (e.g.  /oo)  with  minor  changes. 

4.1  Stability  and  /2-Stability 

Definition  2  Let  (W,E,G)  be  an  implicit  linear  system  over  I2  The  system  is  l2-stable  if  G 
is  left  invertible  in  C{l2),  i-e.  3L  G  C{l2)  such  that  LG  =  7. 

An  equivalent  characterization  is 

Proposition  3  Gw  =  0  is  l2-stable  if  and  only  if 

inf{||6'u;||  :  w  £  I2,  ||w||  =  1}  >  0  (16) 


Interpreting  the  definition,  /2-stability  implies  that  the  I2  behavior  B  of  the  system  is  the 
trivial  space,  and  that  this  property  is  not  “sensitive”  to  equation  error:  an  arbitrary  small 
equation  error  e  does  not  allow  solutions  w,  ||u;||  =  1,  to  the  equation  Gw  =  e. 

An  important  remark  is  that  in  this  framework,  stability  is  only  a  property  of  systems 
with  no  free  variables  (“autonomous”  systems,  in  the  language  of  [19]),  i.e.  with  at  least  as 
many  effective  equations  as  variables  in  w,  so  that  only  the  trivial  behavior  is  left  inside  I2. 

We  will  now  compare  this  definition  with  the  standard  one,  by  considering  the  feedback 
interconnection  of  Figure  3.  The  maps  M,  N  are  (possibly  non-causal)  operators  on  I2. 

This  interconnection  can  be  represented  by  the  following  equations  in  the  variables  m, 


— “1 

1 

_ 1 

w 

i 

_ 1 

z 

(17) 
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According  to  Definition  2,  /2-stability  implies  left  invertibility  of  G 


I  -M 
-N  I 

is  slightly  weaker  than  the  usual  definition,  which  requires  G  invertible  {G  left  invert! 


.  This 
ale  and 

also  G  onto  E).  In  other  words,  Definition  2  is  weaker  in  the  fact  that  equation  errors  are 
not  required  to  be  free  to  vary  over 

The  reason  for  this  weakened  definition  in  the  case  of  implicit  systems  is  that  we  want 
to  extend  the  notion  to  systems  which  are  over-constrained  (more  equations  than  variables), 
such  as  the  example  considered  in  (13-15).  In  these  cases,  the  operator  will  not  be  onto  in 
general,  and  this  should  not  be  required:  the  equation  errors  need  not  be  “free”  since  they  are 
not  physical  noises  (which  should  be  included  in  the  w  variables);  they  just  provide  a  means 
of  testing  sensitivity  of  equations. 

We  now  turn  to  the  notion  of  stability  for  causal  systems.  We  consider  the  following 
definition,  which  extends  the  standard  theory  [22]: 


Definition  3  Let  (W,E,G)  be  an  implicit  linear  system  over  l2e:  W  =  lie,  ^  =  /fe,  G  : 
W-^E  causal  map  .  The  system  is  said  to  be  well  posed  if  G  has  a  causal  left  inverse 
L  :  E^W.  The  system  is  stable  if  it  is  well  posed  and  in  addition,  L  has  finite  gain  on  the 
range  of  G,  i.e.  Bq  G  R"*"  such  that  for  every  e  =  Gw,  and  every  T,  ||Pt^'^||  < 


Remark:  An  immediate  consequence  of  the  definition  is  that  a  stable  system,  with  errors 
e  G  h,  only  allows  signals  w  £  I2  to  satisfy  the  equation  Gw  =  e. 

Definition  3  can  also  be  interpreted  in  terms  of  the  infinite  matrix  representation  (46)  of 
the  equation  operator  G,  and  its  truncations  G^  as  in  (2). 

Proposition  4  The  implicit  system  (W,E,G)  over  l2e  is  stable  if  and  only  if 

info:(6'^)>0  (18) 


To  illustrate  these  definitions,  we  will  consider  the  case  of  an  “autoregressive”  system  (see 
[19])  defined  by  linear  time-invariant  equations  of  the  form  G  =  R{X),  where  R{^)  is  a  p  X  g 
polynomial  matrix. 
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Proposition  5  Let  be  a  p  X  q  polynomial  matrix.  Then 

(i)  The  system  R(X)w  =  0  is  l2-stable  if  and  only  if  R(()  has  full  column  rank  for  all  ^  in 
the  unit  circle  =  1) 

(ii)  The  system  R(X)w  =  0  is  stable  if  and  only  if  R{^)  has  full  column  rank  for  all  ^  in 
the  unit  disk  <  1) 

Remark:  The  same  argument  carries  through  if  R{^)  is  a  rational,  rather  than  polynomial 
matrix,  with  no  poles  on  the  unit  circle  (respectively  the  unit  disk). 

Example  1  Consider  G{X)  =  1  —  2A.  Setting  W  =  E  =  l2e,  we  find  that  the  implicit  system 
over  l2e  is  not  stable,  since  the  signal  w{t)  =  2*,t  >  0,  which  is  in  /2e\^2  gives  e  =  Gw  G  h 
(e(0)  =  1,  e{t)  =  0  for  t  >  0). 

However,  the  system  over  I2  (with  W  =  E  =  I2)  is  l2-stable  since  inf|[i„||=i  ||G'tn||  =  1  >  0. 

4.2  Robust  Stability  and  /2-Stability 

In  this  section  we  consider  the  case  of  implicit  uncertain  systems  as  in  (4).  For  simplicity  we 
will  assume  henceforth  that  the  nominal  equation  map  M  is  in  C{l2)  for  the  /2-stability  case, 
or  in  Cc{he)  for  the  stability  case. 

Definition  4  Let  M  G  Eih)}  ^  C  >C(/2)-  The  implicit  system  (4)  has  robust  l2-stability  if  it 
is  l2-stable  for  each  A  G  Ba- 

Robust  /2-stability  implies  that  for  each  A  G  Ba,  ¥’(A,  M)  has  a  bounded  left  inverse.  In 
some  treatments  (e.g.  [17])  a  uniform  notion  of  robust  stability  is  employed,  which  in  this 
context  would  imply  that  there  is  a  uniform  bound  on  the  norms  of  the  left  inverses  across 
B A ;  we  will  not  pursue  this  refinement  here  since  it  is  appears  to  be  no  stronger  in  most  cases 
of  interest 

The  constrained  I/O  robust  performance  problem  posed  in  Section  3.4  converts  to  robust 
/2-stability  of  the  corresponding  implicit  system:  a  negative  answer  to  Q  implies  a  trivial 
^We  acknowledge  G.  Dullerud  for  pointing  this  out;  see  also  the  proof  of  theorem  10. 
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/2-kernel  for  every  A.  Robust  I2  stability  is  slightly  stronger  in  requiring  this  property  to  be 
insensitive  to  equation  errors,  but  this  is  desirable  anyway.  This  example  will  be  reconsidered 
in  5.3. 

The  following  proposition  provides  a  simplified  representation: 


Proposition  6  The  implicit  system  (4)  has  robust  l2-stabiUty  if  and  only  if 

inverse  L. 

z  =  Q,  where  A  =  A  —  BLC,  C  —  C  —  DLC,  has  robust 


(i)  D  is  I2- stable,  with  bounded  left  inverse  L. 

I  —  AA 

(a)  The  implicit  system  ' 


I2 -stability. 


C 


The  previous  result  has  reduced  robust  I2  stability  of  (4)  to  a  nominal  /2-stability  condition 
(i)  plus  a  robust  /2-stability  condition  (//)  in  a  simplified  setup.  We  now  turn  to  the  situation 
of  robust  stability  in  systems  over  l2e- 


Definition  5  Let  M  6  £c(4e);  A  G  Tc(4e)-  The  implicit  system  (4)  has  robust  stability  if  it 
is  stable  for  each  A  G  Ba  • 


As  is  natural  from  the  definition  of  stability,  robust  stability  means  that  for  each  A  G  Ba, 
(p(A,M)  admits  a  causal  left  inverse  which  has  finite  gain  on  the  range  of  (p(A,M).  For 
the  reduction  result  corresponding  to  Proposition  6,  the  nominal  stability  condition  will  be 
slightly  strengthened  to 

(/')  D  has  causal  left  inverse  X,  which  has  finite  gain  over  l2e- 

The  strengthening  comes  from  the  fact  that  the  finite  gain  of  L  is  required  over  all  /2e,  not 
just  the  range  of  D.  This  does  not  appear  to  be  a  major  limitation,  since  it  is  no  stronger  for 
the  case  when  D  is  LTI,  as  shown  in  the  proof  of  Proposition  5. 

Proposition  7  Consider  the  system  (4)  where  A 
Then  the  system  is  robustly  stable  if  and  only  if 
i  =  A  -  BLC,  C  =  C  -  DLC. 


G  T-dhe),  and  M  G  Ccihe)  satisfies  (i'). 

l-aA^ 


c 


0  is  robustly  stable,  where 
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The  previous  results  show  that  for  robustness  analysis,  it  suffices  to  consider  the  canonical 
implicit  system 

where  for  simplicity  we  have  renamed  A^C  as  A,C.  This  case  will  be  considered  in  the  rest 
of  the  paper. 

It  is  useful  to  compare  this  setup  with  the  question  of  robust  stability  in  standard  robust 
control,  which  specifies  the  invertibility  of  (/  —  AA).  The  main  difference  is  that  (19)  allows 
for  additional  constraints  defined  by  the  C  equations.  A  problem  with  more  equations  than 
variables  such  as  the  one  considered  in  (13-15)  will  result  in  the  presence  of  these  additional 
equations. 

A  problem  where  the  C  equations  do  not  appear  will  be  termed  the  “unconstrained”  case. 
For  example,  if  the  IQCs  (14)  are  eliminated  from  (13-15),  the  problem  can  be  reduced  to  the 
equations 

which  are  in  the  standard  form  (/  —  AA)z  =  0. 

In  an  unconstrained  problem,  the  only  apparent  difference  with  the  standard  case  is  the 
fact  that  our  definition  of  stability  only  specifies  left  invertibility  of  (/  -  AA).  The  following 
proposition  shows  that  the  two  notions  become  in  fact  equivalent  when  they  are  considered 
across  Ba. 

Proposition  8  Let  A  G  £(^2),  A  C  Then  (7  — AA)  has  bounded  left  inverse  VA  G  Ba 

if  and  only  if  (I  —  A  A)  has  bounded  inverse  VA  G  Ba- 

A  similar  result  can  be  obtained  for  the  case  of  robust  stability  of  systems  over  l2e-  This 
shows  nothing  changes  when  using  the  left-invertibility  notion  for  the  unconstrained  case;  on 
the  other  hand,  this  is  the  only  reasonable  notion  for  constrained  problems. 
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5  Robustness  Analysis  with  LTV  uncertainty 

In  this  section  we  will  focus  on  the  class  of  implicit  uncertain  systems  as  in  (19)  where  equation 
maps  A,  C  are  linear  time  invariant  (LTI)  and  the  uncertainty  set  A  consists  of  arbitrary 
linear- time  varying  perturbations  with  spatial  structure  as  in  (3). 

The  main  result  in  this  section  is  a  necessary  and  sufficient  condition  for  robust  stability. 
This  result  extends  the  “scaled  small  gain”  sufficient  conditions  for  robust  stability  in  the 
standard  input  output  framework,  and  recent  results  on  the  necessity  of  the  constant  scales 
tests  for  linear  time-varying  perturbations,  obtained  by  Shamma  [17]  and  Megretski  [9,  10] 
(and  also  previously  in  the  setting  by  Khammash  [8]). 

Once  more  the  I2  and  l2e  settings  are  considered  separately;  we  will  concentrate,  however, 
in  the  I2  setting  since  the  main  motivation  for  over-constrained  problems,  given  in  section  3.4, 
refers  to  this  case.  These  issues  will  be  discussed  further  in  section  5.3. 

5.1  A  Necessary  and  Sufficient  Condition  for  Robust  /2-Stability 

We  now  consider  analysis  in  the  I2  setting;  A  ,C,  A  are  in  £(/2)  of  appropriate  dimensions 
(over  Z  or  Z+,  and  are  not  restricted  to  be  causal);  A  ,C  are  LTI. 

In  view  of  the  duality  between  IQCs  and  implicit  LFT  descriptions  which  was  shown 
in  section  3.3,  the  theorem  given  below  is  a  version  in  the  implicit  LFT  framework  of  the 
S-procedure  losslessness  results  of  Megretski  and  Treil  [10].  The  main  extension  needed  is 
to  capture  the  61  blocks,  or  equivalently,  to  consider  matrix-valued  IQCs;  this  is  done  by 
extending  the  “V  set”  method  in  standard  /i- analysis  [11].  This  extension  has  more  than 
purely  theoretical  interest  since  these  representations  arise  naturally  in  the  context  of  set 
characterizations  of  white  noise,  as  shown  in  Part  11. 

We  begin  with  some  further  notation.  For  a  delta  structure  A  of  the  form  (3),  let  Y  be 
the  set  of  constant,  hermitian  scaling  matrices  Y  that  commute  with  the  elements  in  A, 

Y  =  {diag  [Yi,. .  .,Yl,  VL+iImi ,  ■  •  • ,  yL+Flmp] ,  Yi  =  Yi*  ,i  =  1 . . .  L,  ul+j  E'R,j=  1 . . .  F}  (21) 
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Y  is  a  real  vector  space,  and  we  can  define  an  inner  product 

(y,y)  =  ^tr(Ky,)  (22) 

8=1  i—\ 

Two  important  subsets  of  Y  are  X  =  {y  €  Y,y  >  0},  and  X  =  {y  G  Y,y  >  0},  the  set  of 
positive  and  nonnegative  scalings,  respectively.  They  are  both  convex  cones  in  Y. 

Given  a  vector  z{t)  G  ^  ^^e  block  structure  introduces  a  natural  partition  of  2, 

Z  =  [zj  .  .  .  ^L+l  •  •  •  %+f] 

Given  an  LTI  map  A  G  C{l2),  an  analogous  notation  is  used  for  the  partition  of  Az  E  I2. 
Consider  the  following  quadratic  functions  of  z  G  ;  where  ^  =  Az: 

00 

t—  —  oo 

00 

=  \\(A^)L+jf  -  \\^L+j\f  =  Y  ^L+jity^L+jit)  -  ZL+j{tyZL+j(t)  j  =  1...F 

tzz—00 

A(z)  =  diag[$i(z),...,$i(z),  0-^+1  (z)/„,, ..  .,o-i+2p(z)/„jp]  (23) 

Now  let  C  G  C(l2)  be  LTI  and  e  >  0.  Define  a  subset  of  Y, 

=  {A(z)  :  ||z||  =  1,  ||Cz||  <  e}  (24) 

Lemma  9  The  closure  ofV  is  convex  and  compact  in  Y. 

Theorem  10  Let  A,C  be  LTI  stable  systems,  as  in  (24).  Assume  A  is  the  set  of  struc¬ 
tured,  otherwise  arbitrary  linear  operators  in  l^.  The  following  are  equivalent: 

(i)  The  implicit  uncertain  system  (19)  has  robust  I2  —  stability. 

(ii)  3  e  >  0  such  that  n  X  =  0 

(m)  3X  €XsuchthatA*XA-X -C*C  <0  (25) 

In  (Hi)  above,  condition  (25)  is  of  the  form  $  <  0  where  4^  is  a  self-adjoint  operator  on 
T]  this  must  be  interpreted  as  a  strong  version  of  negative  definiteness,  u)  <  — p||m||^ 
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for  some  p  >  0.  For  the  case  of  A,C  finite  dimensional,  this  condition  can  be  tested  in  the 
frequency  domain  in  the  form 

3  Z  G  X,  such  that  A(e^‘")*XA(e'“)  -  X  -  C(e^‘^yC{e^^)  <  0  Vw  G  tt]  (26) 

This  test  is  a  Linear  Matrix  Inequality  (LMI,  [1])  evaluated  over  frequency,  and  lends  itself 
to  available  convex  optimization  tools. 

5.2  Necessary,  Sufficient  Conditions  in  the  /2e  setting 

As  remarked  before,  the  motivation  we  have  provided  for  constrained  problems  refers  mainly 
to  the  I2  setting.  It  is  well  known,  however,  that  in  the  standard  unconstrained  case,  the 
condition  A*XA  —  X  <  0,  (or  equivalently  XiAX~i  <  1)  is  sufficient  for  robust  stability 
in  the  l2e  sense  (this  is  a  consequence  of  the  small  gain  theorem  [22]).  Also,  [17,  9]  show  it  is 
a  necessary  condition. 

It  therefore  seems  natural  to  explore  this  issue  in  the  constrained  case  of  (19).  We  first 
consider  necessity: 

Theorem  11  Let  A,  C  be  LTI  in  jC-dhe),  ^  is  the  set  of  structured,  otherwise  arbitrary 
causal  linear  operators  in  jCdhe)-  If  (19)  over  l2e  is  robustly  stable  then  (25)  holds. 

For  the  proof,  it  clearly  suffices  to  prove  condition  (u)  in  Theorem  10;  the  only  required 
modification  in  the  argument  presented  for  Theorem  10,  is  that  causality  of  the  destabilizing 
perturbation  must  be  ensured.  This  can  be  done  following  the  lines  of  the  proof  in  [17]  for 
the  I/O  case;  for  reasons  of  brevity  this  construction  will  not  be  developed  here. 

Regarding  the  issue  of  sufficiency  of  condition  (25),  we  consider  the  following  example. 

Example  2  Let  A  =  2,  C  =  1  —  2A.  A  G  Cdhe)-  For  any  X  G  R"*", 

{{A*XA  -X-  C*C)v,v}  =  3X  Ijujf  -  ||(1  -  2X)v\f  <  (3X  -  1)  ||u||^  (27) 

so  for  0  <  X  <  I  (25)  is  satisfied.  But  the  perturbation  A  =  A  gives 
which  is  unstable,  as  shown  in  Example  1. 
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This  example  shows  that  in  general,  condition  (25),  (or  (26))  is  not  sufficient  for  robust 
stability,  even  for  time  invariant  perturbations:  (25)  does  not  provide  information  of  the 
behavior  outside  I2.  One  could  think  of  strengthening  condition  (26)  to  include  frequency 
points  inside  the  disk;  this  would  eliminate  the  previous  counterexample,  and  in  fact  guarantee 
robust  stability  with  LTI  perturbations  (see  Part  II),  but  it  is  still  not  sufficient  for  the  LTV 
case.  A  counterexample  for  this  is  the  system  A  =  4A  —  8A^,  C  =  1  —  4A  +  4A^;  we  omit  the 
rather  lengthy  verification 

We  can  state,  however,  a  partial  result  which  is  applicable  in  many  cases. 

Theorem  12  Let  A  be  LTI  in  £c(^2e)j  and  C  be  a  static  map.  If  (25)  holds,  then  system  (19) 
is  robustly  stable  over  the  class  of  structured,  otherwise  arbitrary  operators  in  Cdhe)- 

Although  this  result  is  not  a  major  extension  of  the  standard  small  gain  theorem  of  [22], 
since  only  static  additional  constraints  are  allowed,  it  is  still  quite  rich  from  the  point  of  view 
of  posing  constrained  robustness  analysis  problems,  as  will  be  shown  below. 

5.3  Analysis  for  Robust  Performance  Revisited 

To  conclude  Part  I  of  this  paper,  we  review  the  example  considered  in  section  3.4  using  the 
results  of  this  section.  Referring  to  (13-15)  and  Figure  2  we  assume  that  H,  P  and  Q  are 

causal,  LTI,  stable  maps.  We  recall  from  section  3.4  that  without  loss  of  generality  in  the 

analysis,  P  can  be  chosen  to  be  a  static  map  kl.  Also  the  perturbations  Ap,  Ac  vary  in  the 
class  of  arbitrary  time- varying  operators. 

The  reduction  procedure  from  Propositions  3,  4,  yields  a  system  of  the  form  (19),  where 
Au  0  0  Hu  Hi2  0 

A  =  0  Ap  0  ,  A  =  H21  H22  0  ,  C  =[0  P  -l],  z=  zp  (28) 

_  0  0  AcJ  [  0  g  oj  L^e. 

If  P  is  static,  we  observe  that  C  is  static,  as  in  the  assumption  of  Theorem  12.  Therefore 
if  A,,  is  LTV  uncertainty,  condition  (25)  is  necessary  and  sufficient  for 

^This  counterexample  can  be  related  to  an  analogous  situation  in  [20],  involving  a  frequency  domain  con¬ 
dition  over  the  right  half  plane. 
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•  robust  /2-stability  in  the  case  A^,  Ap,  Ac  in  £(/2), 

•  robust  stability  in  the  case  A^,  Ap,  Ac  in  Ccihe), 

If  we  review  the  problem  statement  that  let  to  equations  (13-15),  the  I2  version  seems  more 
appropriate,  but  we  might  also  be  interested  in  a  “hybrid”  problem,  by  testing  strong  stability 
in  the  variables.  More  precisely,  A„  may  be  considered  a  causal  perturbation,  and  it  must 
be  guaranteed  in  the  first  place  that  the  variable  does  not  “blow  up”  (it  remains  in  I2  if 
d  and  the  equation  errors  are  in  I2).  Once  this  is  known,  the  analysis  can  be  restricted  to  I2 
and  the  (possibly  non-causal)  perturbations  Ap,  Ac  can  be  considered,  casting  the  robust 
performance  analysis  as  a  robust  I2  stability  question,  as  argued  in  section  3.4. 

This  hybrid  question  is  also  answered  by  the  test  (25);  in  fact,  since  the  first  block  of  C  is 
zero,  the  upper  portion  of  A* X A  —  X  —  C*C  is  Hl^X^Hn  — which  provides  the  standard 
robust  stability  test  in  the  z^,  variables,  in  addition  to  the  robust  I2  stability  test  on  all  the 
variables. 

6  Conclusions 

The  work  reported  in  Part  I  of  this  paper  provides  the  foundation  for  a  more  general  robustness 
analysis  theory,  which  extends  the  standard  theory  based  on  the  small  gain  theorem.  In  this 
approach,  we  abandon  the  concepts  of  “input-output  maps”  and  “gains”  in  favor  of  equations 
and  signal  constraints^  and  the  central  analysis  question  is  to  test  whether  there  exist  solutions 
to  these  equations. 

The  results  in  this  paper  demonstrate  that  nothing  is  lost,  from  a  mathematical  point  of 
view,  by  adopting  this  approach  for  ana/j/s/s  instead  of  the  standard  input-output  formulation; 
on  the  contrary,  the  analysis  setup  presented  in  section  3.4,  further  developed  in  Part  II,  shows 
evidence  of  substantial  advantages. 

There  are  still  reasons  to  preserve  the  standard  “signal-flow”  approach,  which  has  led  to 
a  large  body  of  knowledge,  since  some  of  its  intuitive  value  for  design  is  lost  in  the  implicit 
formulation.  The  conclusion  is,  however,  that  if  research  is  not  confined  to  the  traditional 
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paradigm,  the  potential  of  the  resulting  theory  will  be  enhanced. 
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Appendix:  Proofs 

Lemmas  1  and  2 

The  only  non-trivial  implication  is  (ii)  (Hi)  in  Lemma  2. 

If  Vi,...Vr,  is  an  orthonormal  basis  of  the  subspace  of  spanned  by  the  coordinates 
Vi, . .  .v^  of  u,  we  write  v  =  Pv,  where  P  is  an  invertible  matrix,  v  =  [hi, . .  0, . .  .0]'  .  Let 

h  =  Pz,  then  (ii)  implies  i  =  [ii, . .  0, . .  .0]',  and 

\\riiVi  +  .  .  .  +  TlrVrW  >  \\riiZi  +  .  .  .  +  Tj^ZrW  V??! ,  .  .  .  7/,.  €  C  (29) 

Now  define  6  :  u  then  ^  :  h,  i,-, i  =  1 . .  .r,  so  div  =  z  which  implies 

8Iv  =  2;.  Also,  by  (29)  and  the  Bessel  inequality,  ||^7i||  <  \\Yl^i=i{u.,Vi)vi\\  <  ||m||.  So  ||h||  <  1. 

□ 

Proposition  3 

For  the  necessity,  note  that  since  LG  =  /,  ||L||  HGrcH  >  ||7n||.  For  the  sufficiency,  condition 
(16)  implies  that  G  is  injective,  and  its  range  is  closed.  An  application  of  the  open  mapping 
theorem  [4]  implies  that  G  has  a  bounded  inverse  on  its  range,  which  can  be  extended  to  a 
bounded  operator  on  E,  resulting  in  a  left  inverse  for  G. 


□ 
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Proposition  4 


[Necessity]  the  definition  of  stability  gives  ^  ||^||  from  which  (18)  follows. 

[Sufficiency]  Let  e  =  infj’>oCT(G^)  >  0.  Then  a  causal  left  inverse  L  can  be  constructed 
by  induction  over  T:  given  a  lower  triangular  left  inverse  of  ,  consider  the  matrix 

'  (~<T  n  1 

^  ^  •  Since  >  0,  Gt+i,t+i  has  left  inverse  Lt+i,t+i,  then 


G  G- 


T+1,T+1 


0 

—Lt+i,t+iGL^  Lt+i,t+i 
is  a  left  inverse  of  Gt+i-  This  procedure  produces  a  causal  operator  X  :  E-^W.  If  e  =  Gw, 
then  Pye  =  G^ PtW  therefore  ||PTe||  >  cUPyXeH  so  X  has  finite  gain  in  the  range  of  G. 


□ 


Proposition  5 

(i)  If  R(e^‘^°)w  =  0,  ro  7^  0,  it  is  easy  to  construct  R  signals  with  norm  1,  in  the  direction 
w,  and  spectrum  supported  in  [o^q  —  |,n;o  —  i],  so  that  ||P(A)rc(*)||  0, violating  /2-stability. 

Conversely  if  R{^)  has  full  column  rank  on  the  unit  circle,  min^^  £(^(6^“))  >  0  and  it  is  easy 
to  construct  pointwise  an  £00  left  inverse  for  R(X). 

(ii)  If  P(0)  has  a  kernel,  then  P(A)  cannot  have  a  causal  left  inverse,  so  well  posedness 
fails.  If  R{^o)w  =  0,  0  7^  l^ol  <  1  then  the  signal  w{t)  =  {~yw,  /  >  0  is  in  l2e\l2,  but 
R{X)w  €  I2,  violating  stability. 

If  R{^)  has  full  column  rank  over  the  unit  disk,  then  the  theory  of  coprime  factorizations 
over  the  stable  ring  (see  e.g.  [18])  implies  that  R(^)  has  a  left  inverse  in  TZHco-  This  is  a 
causal  left  inverse  which  has  finite  gain  (over  all  l^e)  so  stability  is  satisfied. 

□ 


Proposition  6 


[Necessity]  If  (4)  is  robust  /2-stable,  setting  A  =  0,  <^(0,  M)  = 


X  0 
C  D 


has  bounded  left 
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inverse,  therefore  D  has  bounded  left  inverse  L.  Fix  A  G  Ba;  the  following  identity  holds: 


(p(A,M)  U  =  V  cp(A,M) 
where  U  and  V  are  the  invertible  operators 
U  = 

and 


(30) 


'  I  o‘ 
-LC  I 

,  v  = 

'  I  -ABL 

0  I 

,  u-^  = 

I  o' 
LC  I 

’  I  ABL 

0  I 

M 


A  0 
C  D 


,  (^(A,M)  = 


/-AA  0 
C  D 


(31) 


(32) 


Also  denote  '9i(A,M) 


.  If  r(A)  is  a  bounded  left  inverse  for  (p(A,M),  then 


I  -  AA 
C 

U~^T{A)V  is  a  left  inverse  for  (p{A,M),  which  implies  $(A,  M)  is  left  invertible. 
[Sufficiency]  If  T(A)  =  Ti  T2  is  a  bounded  left  inverse  for  W(A,M)  then  (note  that  by 
definition  of  C,  LC  =  0), 


.  1 

to 

i-^ 

1 

b 

' l-aA  o' 

- > 

0 

_ 1 

C  D 

I 


(33) 


therefore  (p(A,M)  is  left  invertible,  and  so  is  ^(A,  M)  by  (30). 


□ 


Proposition  7: 

Defining  U,  V  as  in  (31),  by  assumption  {i')  P,  are  in  Cdhe)-  From  (30), 


robust  stability  of  (4)  reduces  to  that  of  (p(A,M)  =  0.  Clearly  this  implies  stability  of 

.  ^  . 

^{A,M)z  =  0,  It  remains  to  show  the  converse  implication. 

Fix  A  G  Ba,  let  [  Ti  T2  ]  be  a  causal  left  inverse  of  $(A,M),  with  finite  gain  on  its 


range,  and  construct  a  causal  left  inverse  for  (p(A,M)  as  in  (33).  The  bottom  portion 
has  finite  gain  by  (*').  The  top  portion  is 


Ti  T2{I-DL)  =  Ti  T: 


I  0 
0  {I-DL) 


0  L 


(34) 
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where 


has  finite  gain  and  maps  a  vector  in  the  range  of  M),  to  the  range 

of  $(A,M)  where  Ti  T2  has  finite  gain,  therefore  the  top  portion  has  finite  gain  in  the 
range  of  y>(A,  M). 

□ 

Proposition  8 

The  result  follows  from  some  results  in  spectral  theory  of  operators.  For  simplicity  we  consider 
the  complex  field  F  =  C,  but  a  similar  proof  can  be  written  for  the  real  case.  For  an  operator 
M  :  /2(C^)-^/2(C®),  the  spectrum  is  defined  as  c(M)  =  {C  G  C,{(I  —  M)  not  invertible}.  A 
subset  of  the  spectrum  is  the  approximate  point  spectrum,  defined  as 

=  {C  G  C,  ((■/  —  M)  not  left  invertible}  =  {C  G  C,  inf  \\((I  —  M){z)\\  —  0}  (35) 

ibil=i 

It  is  known  (see  [4])  that  the  boundary  of  the  spectrum  is  in  the  approximate  point  spectrum. 
Therefore  the  spectral  radius  p{M)  =  max{|C|  :  (  G  f(M)}  is  achieved  at  Co  G  ‘;ap{M). 

We  now  prove  the  only  if  portion  of  the  proposition  (the  other  direction  is  trivial).  Assume 
there  exists  Aq  G  Ba,  such  that  /  -  AqA  not  invertible.  Therefore  p(AoA)  >  1.  This  implies 
there  exists  Co  G  ?ap(AoA),  with  |Co|  >  1-  Therefore  Co^f  —  AqA  is  not  left  invertible,  which 
gives  I  —  AA  not  left  invertible,  with  A  =  ^A  G  Ba,  a  contradiction. 

□ 

Lemma  9 

Since  A  G  C{l2),  V'  is  bounded,  therefore  is  compact.  For  the  convexity  proof,  it  suffices 
to  show  that  co(V')  C  V%  where  co(V^)  is  the  convex  hull  of  VG  Consider  the  convex 
combination  ah.{z)  +  (1  —  Q!)A(t;)  of  two  elements  of  (a  finite  number  of  terms  can  be 
handled  in  a  similar  way);  (||^||  =  ||n||  =  1,  IIC^H  <  e,  \\Cv\\  <  e,  0  <  a  <  1).  Define 
yik)  _  (^a)^z  +  (1  —  a)^X'°v.  Then 

=  a\\z\f  {1  —  a)\\v\^  —  a)^ {z,X^v)  (36) 


I  0 
0  (I~DL) 
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=  a\\Czf +  {l-a)\\Cv\f +  2{aY^{l-ay^{Cz,\^Cv)  (37) 

where  (37)  uses  the  time  invariance  of  C.  The  last  terms  in  (36),  (37)  have  limit  0  as  k-^-oo. 
Therefore  a  +  (!-«)  ||u||^  =  1,  and  ||C’u(*'^|l  <  e  for  large  k,  and  therefore 

C ||"(fc)||  <  e  foT  k  >  ko.  A  similar  argument  using  the  time  invariance  of  A,  shows  that 

$,(y(*))  a^i(z)  +  (1  -  o)$i(u)  i  =  l...L  (38) 

aL+j{v^'‘^)  aai+jiz)  +  {1  -  a)aL+j{v)  j  =  l...F  (39) 

This  gives  A(||^)  aA(z)  +  (1  -  q;)A(u)  G 

□ 

Theorem  10 

(i)  (ii)  By  contradiction,  assume  that  for  all  e  >  0,  fi  X  ^  0.  For  a  fixed  e  >  0,  we  can 
therefore  find  z  E  I2  such  that  ||^||  =  1,  \\Cz\\  <  e,  and  A(^)  +  e^I  >  0.  Let  ^  =  Az.  We  have 

||a+i|p-|kL+il|^  +  c^>0  j  =  l...F  (40) 

00 

>  0  i  =  l...L  (41) 

tz=--00 

Focusing  on  (40),  a  slight  extension  of  Lemma  1  shows  that  there  exists  an  LTV  operator 
Ai+j,  ||Ai+j||  <  1  and  an  error  signal  e^+j,  ||ei+j||  =  0(e)  such  that  AL+jU+j  +  ^L+j  =  zl+j. 

A  similar  argument  (extending  Lemma  2)  is  used  for  (41).  The  result  is  a  structured  LTV 
operator  A,  ||A||  <  1  and  an  error  signal  e,  ||e||  =  0(e)  such  that  (I  —  AA)^  =  e.  Therefore 

z  =0(e).  (42) 

Since  this  holds  for  any  e,  there  exist  sequences  of  signals  z^  and  perturbations  A^,  such  that 

||^'^||  =  1,  ^-^0  (43) 

It  remains  to  show  that  a  single  perturbation  A  and  a  sequence  z^  can  be  found  with  the 
same  property.  The  following  construction,  (which  shows  that  the  concept  of  uniform  robust 
stability  is  no  stronger  in  this  case)  is  due  to  Dullerud  [7].  For  simplicity  we  consider  T  =  Z"*". 
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Starting  from  (43),  let  ||(/—  Aj;A)z*||  <  e*  with  e*,  0.  For  fixed  k  we  can  find  a 

truncation  time  Tj.  such  that 

11(1-  <  €,  and  \\{I  -  P'^>’)Az’^\\  <  (44) 

This  follows  from  the  fact  that  P'^Af.P^v  for  any  v;  P^A^P'^  is  identified  with  the 

truncated  matrix  A^  as  in  (2).  Now  define  the  operator  A  with  infinite  matrix  representation 

'Al”  0  0  •••■ 

0  Af^  0  •  •  • 

A  =  :  ;  ;  •  •  •  (45) 

0  Al’’ 

Let  Tk  be  defined  by  Tq  =  0,  r^+i  =  r*  +  T*  for  k  >  0,  and  z’^  =  .  From  the  time 

invariance  of  A,  (I  —  AA)z’‘  =  -  AA’’’'^^*  =  A’‘‘(/  -  A’‘A)z’'  where  A*  has  infinite 

matrix  representation 

A^'‘  0  0 

0  (46) 

It  follows  that 

\\{I  -  AA)z’‘\\  =  (/-  P^'‘A;tP^M)z*  +  (P^'‘AfcP^‘  -  A*)A^H 

<  11(7-  P'^^AtP^M)^"^!  +  1(1-  P^'')^^*^  II  <  2e,  0  (47) 

is  not  /2-stable. 

□ 

(n)  =^>  (Hi) 

and  X  are  disjoint  convex  sets  in  the  inner  product  space  Y,  is  compact  and  X  is 
closed.  We  can  use  a  hyperplane  separation  argument  to  find  X  G  Y,  77  >  0  such  that 

(X,A)<a-77<a<  (X,y)  VAe  V%  YgX  (48) 


Since  ||;2r*||  =  1,  and  ||C2r^ 


0  from  the  time  invariance  of  C, 
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Since  X  is  a  cone,  a  can  be  chosen  to  be  0.  It  is  easy  to  show  that 

(X,  y)  >  0  VF  G  Y,  y  >  0  x  >  o  (49) 

A  small  perturbation  of  X  ensures  X  >  0  (therefore  X  G  X),  and  by  continuity  and  com¬ 

pactness  of  we  can  modify  77  to  achieve 

(X,A)< -7/<  0  VA  G  (50) 

for  the  new  X.  Furthermore,  by  scaling  down  X  (and  rj)  we  can  ensure  (50)  holds  and 

||A*XA-X||  =  7"<6".  (51) 

Now  for  any  2;  =  I2,  ||-^||  =  1,  ^  =  Az, 


L  F 

(X,A(z))  =  ^tr(Xi$i(z))  + ^X£+j(7i+j(z) 

i=l 

=  {XAz,Az)- {Xz,z)  =  {{A*XA- X)z,z) 


(52) 


Let  ^  =  A*XA-  X  -  C*C. 

If  lie'll!  <  e,  then  A{z)  G  V'^  so  {^z,  z)  <  {{A*XA  —  X)z,  z)  =  (X,  A(2;))  <  -7  <  0  by  (50). 
If  ||C2:||  >  e,  then  z)  <  ((A*XA  —  X)z^  z)  —  HCzH^  <  7^  —  <  0  using  (51). 

Therefore  z)  <  max{—r],'y^  —  e^)  <  0  for  ||z||  =  1,  which  implies  $  <  0. 

(in)  =7  (i) 

Fix  X  >  0  which  solves  (25).  Without  loss  of  generality  it  can  be  assumed  that  X  =  /. 
This  results  from  the  fact  that  invertible  operations  yield 


X5  0 
0  / 


I- AA 
C 


X^ 


/-  AA 
C 


(53) 


where  A  =  X ^AX~^,  C  =  CX  2  verify  (25)  with  X  =  J.  The  notion  of  negative  definitenes 


allows  by  continuity  to  find  <  1  such  that 


A*A-a2/-C*C'<  0. 


Let  e,z  e  I2  (ll^^ll  =  1)  and  A  (||A||  <  1)  satisfy 


I-  AA 
C 


z  =  e  = 


ei 

62 


(54) 


(55) 
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(56) 


IIAA^IP  <  llA^ll'  =  {A*Az,  z)  <  a\z,  z)  +  {C*Cz,  2)  =  +  Heslf 

Therefore  ||AA;2||  <  a  +  ||e2||,  and  ||ei||  =  \\z  —  AA2;||  >  (1  —  a)  —  ||e2||  This  shows  that 
IHi>(l|eill  +  l|e2ll)/^>(l  -  q.)I'/2,  which  proves  robust  /2-stability  by  Proposition  3. 

□ 


Theorem  12 


By  the  same  argument  given  above,  assume  without  loss  of  generality  that  X  =  I,  and  let  a 
satisfy  (54).  Let  T  be  fixed,  z,e  satisfy  (55),  with  truncation  P'^ z  of  norm  1.  From  causality 


I- AA 
C 


P^z  =  P^e 


P'^ei 

P^e2 


(57) 


Since  P'^  z  G  h,  from  (54)  we  obtain 

\\AP^zf<a^  Wp-^zf  +  WCP^zf 

Since  C  is  static,  CP'^z  =  P^e2;  also  \\P'^ AAP'^ z\f  <  \\AP'^z\\^,  leading  to 

llP^eill  =  l|P^(/-AT)P^z||  >  1-a-  ||P^e2|| 


(58) 


(59) 


We  conclude  that  ||P^e||  >  c  for  a  fixed  c  >  0,  independent  of  T,  A,  which  implies  that 

1 TN 

infr  all  ^  I  I  >  0  for  every  A  giving  robust  stability  from  Proposition  4. 


C 


□ 
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